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Abstract. We derive error estimates for multinomial approximations of Amer- 
ican options in a multidimensional jump-diffusion Merton's model. We assume 
that the payoffs are Markovian and satisfy Lipschitz type conditions. Error 
estimates for such type of approximations were not obtained before. Our main 
tool is the strong approximations theorems for i.i.d. random vectors which were 
obtained in |14| . For the multidimensional Black-Scholes model our results can 
be extended also to a general path dependent payoffs which satisfy Lipschitz 
type conditions. For the case of multinomial approximations of American op- 
tions for the Black— Scholes model our estimates are a significant improvement 
of those which were obtained in [8] (for game options in a more general setup) . 



1. Introduction 

This paper deals with multinomial approximations of American options arbitrage- 
free prices in the multidimensional jump-diffusion Merton's model with finite hori- 
zon. The Merton model is a generalization of the Black-Scholes model and it 
allows the stock to have jumps of compound Poisson type. We consider a Markov- 
ian payoffs which satisfy Lipschitz type conditions and derive error estimates for an 
appropriate multinomial approximations. In the multidimensional Black-Scholes 
model these error estimates can be also derived for general path dependent payoffs 
which satisfy Lipschitz type conditions. 

For American options in finite horizon Merton's model arbitrage-free prices can 
not be calculated explicitly. Since multinomial models are active on a discrete set 
of times and defined on a discrete probability space, then arbitrage-free prices in 
these models can be calculated efficiently by dynamical programming algorithm. 
Thus convergence results for multinomial approximations provides an efficient tool 
to evaluate arbitrage free prices in the Merton model. 

Several papers dealt with multinomial approximations of American options in 
models with jumps (see, [10] and [E])- In both of the papers the authors considered 
the case with one risky asset and used the weak convergence approach. The main 
tool that was used in the above papers is the stability results for Snell's envelopes 
under weak convergence which were obtained in [11) . In [5] the weak convergence 
approach was used to show convergence results for multinomial approximations 
of game options (which were introduced in [6]) in the Merton model. The main 
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disadvantage of the weak convergence approach is that this machinery can not 
provide, in principle, speed of convergence estimates. Thus no error estimates were 
obtained for multinomial approximations of American options in the Merton model. 

Clearly, from practical point of view it is valuable to find estimates of the cor- 
responding errors. In order to obtain error estimates we should consider all the 
market models on the same probability space, and so methods based on strong 
approximation theorems come into picture. Strong approximation theorem allows 
to construct a probability space which contain all the markets models such that 
the risky assets in the discrete time models will be "close" with respect to the sup 
norm to the risky assets in the continuous model. Several authors applied strong 
invariance principles in order to obtain error estimates for American and game op- 
tions in the one dimensional BS model (see, [5], [7] and [3]). In all of these papers 
the authors used the Skorokhod embedding tool of i.i.d. random variables into the 
one dimensional Brownian motion. This tool can not be applied for the multidi- 
mensional case. In [8] the author studied discrete time approximations of Dynkin's 
game values for the multidimensional Brownian motion. The main tool that was 
used there is strong approximation theorems which were developed in [Ij and they 
work for sequences of random vectors with close characteristic functions. 

In 1^14] the author considered a new approach to strong approximations in the 
multidimensional case. He showed that for a given sequence of i.i.d. random vec- 
tors X^^\X^'^K..., and a random vector Y which has the same expectation and 
covariance matrix as X^^\ it is possible to construct a sequence of i.i.d. vectors 
y(i)^y(2) _ such that Y^^^ ^ Y and the normalized sums of the last sequence will 
be "close" to the normalized sums of the first sequence. Furthermore for any k < m, 
Y^'^'> is independent of X^™'. The above approach will be the main tool that we use 
in order to establish the results in this paper. For Merton's model strong approxi- 
mation theorems were not used before. An interesting question which is still open, 
is whether the method from [lAi can be applied for game options approximations. 

Let us remark, that the estimates for the Brownian motion by means of normal- 
ized sums of independent random vectors obtained in jl4j are much better than 
those which were obtained in J^ . 

Main results of this paper are formulated in the next section where we also 
introduce the notations that will be used. In Section 3 we derive auxiliary lemmas 
that we use. In Section 4 we complete the proof of main results of the paper. 



Consider a complete probability space (fi, J^, P) together with a standard d- 
dimensional continuous in time Brownian motion {W{t) = {Wi{t), ...,Wd{t))}'^Q, 
a Poisson process {N{t)}'^Q with intensity A and independent of W, and a sequence 
of i.i.d. random vectors {t/^*^ = {u[^\ C/^^^)}5^i with values in (—1, oo)'', inde- 



E denotes the expectation with respect to P. A d-dimcnsional Merton's model 
with horizon T < oo consists of a savings account which given by 



2. Preliminaries and main results 




(2.1) 



b{t) = 6(0)exp(ri), 6(0), r > 0, < t < T 
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and of d risky stocks {S{t) — {Si{t), Sdit))}J^Q given by 
(2.2) 

d d N(t) 

S,{t) = SM exp((r + - E 4/2)^ + E IT + f^^'^)' 1 < * < ^ 

i=i i=i i=i 

where cr = i'^ij)i<i,3<d is a nonsingular matrix, Si{0) > and without loss of 
generality we assume that 

(2.3) fi, = -\Eu'l^\ l<i<d. 
Consider an American option with the payoff process 

(2.4) Y{t)^F{S{t),t), 0<t<T 

where F : M."^ x ]R+ is a function such that for some constant L > 1 and for 

any t > s > and v,v E M"^, 

d 

>. 

Let T be the set of stopping times with respect to the natural filtration generated by 
S (which satisfies the usual conditions) with values not exceeding T. The equality 
(j2.3[) guaranties that the probability measure P is a martingale measure. Thus the 
term 

(2.6) V = sup E{exp{-rT)Y{T)) 

reT 

gives an arbitrage-free price for the American option. 

For any nG N define a sequence of i.i.d. random vectors {[/"'* = (J7"'', .--i '*)} =i 

by 

M{n) 

(2.7) C/- = Y: i^"^' - l)I{-„-vs_i<^<.)<.„-v._ij, l<J<d,^en 

fc=i ' ^ 

where Iq = 1 if an event Q occurs and = if not, and M(n) e N satisfies 

(2-8) E^(t^f'l{ar)>HM„-r/a_ij) < 



(2.5) \F{v, t) - F{v, s)\<lJ2 W^ - ^.1 + Ht - + E 



Notice that 

(2.9) E\U^''> - < e N, I < j < d. 

Next we describe the discrete time markets which we use to approximate the Merton 
model. Let A e Md+i(R) be an orthogonal matrix such that it last column equals to 
(--^==, --^==). Let = {1, 2, d + 1}°° be the space of infinite sequences uj = 

(wi, W2, •■•); G {1; 2, d+1} with the product probability — {j^i dTi^°° ' 
Define a sequence of i.i.d. random vectors C*-^-*, ... by 



(2.10) ^«(c^) = VdTT{A^a,A^a-,A^.d), i e N. 
Observe that 

(2.11) = and i^Cf ^ 1 if i = j and =0 otherwise. 
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The probability space {ft^, P^) was introduced in For any n we extend (ilj, P^) 
to a probability space {^n,Gn, Pn) such that it contains a three independent se- 
quences of i.i.d. random vectors and {u^''^}^^i- The second 
sequence is a sequence of Bernoulli random variables such that Pn{p^'^ = 1} = 
1 — exp(— AT/rt) and the second sequence satisfies m"'^ ~ C/"'^. 

For any n e N, < fc < n and 1 < i < d set A^"'*^ = P"'™ and 

(2.12) 

i\fc 



= ^,(0)exp(rfcT/n) (l+^f E^'.^j 



j=i (1 + (1 - exp(-Ar/n))£;„<'^) 

where ^ <t < ^ ' and denotes the expectation with respect to P„. Con- 
sider a multinomial n-step market which is active in the moments 0, ^, ...,T 
and consists of a savings account which given by (|2.1I) and of d risky stocks 5*^"^ = 
(Si^\ s'^^^) given by (|2.12p . Next, we introduce an American option with the 
payoff process 

(2.13) y(")(fc) = F{S^''-\kT/n),kT/n), < k < n. 

Let Tn be the set of stopping times with respect to the filtration {a{slp\ s'^^ 5*12: )}fc=o 
with values in {0, 1, ...,n}. Observe that for any n, Pn is a martingale measure for 
the n-step market. Thus 

(2.14) K= sup £;„(exp(-Tr/n)r(")(r)) 

rer„ 

is an arbitrage-free price of the n-step market. The following theorem says that 
the arbitrage-free prices of the n-step markets converge to the arbitrage-free price 
of the Morton model and provides an estimates on the error terms. 

Theorem 2.1. For any e > there exists a constant such that for any n 

(2.15) \V^Vn\<C,n'-^. 



Remark 2.2. Theorem \2.1\ can he extended to a case where we have a finite number 
of Poisson clocks. Namely, consider a complete probability space (H, T, P) together 
with a standard d-dimensional continuous in time Brownian motion W , m inde- 
pendent Poisson processes N^^\ N^"^\ which are independent ofW, and for any 
1 < k < m a sequence of i.i.d. random vectors {U^''^}°°^i with values in (— IjCxd)''. 
We assume that the sequences are independent of each other and independent of W 
and N'-^\...,N'-"'\ We also assume Ep^'^l'^ < oo, for j < d and k < m. The 
risky assets are given by 

d d m N^'Ut) 

(2.16) = SM exp((r -f A^. - ^ 4/2)* + E '^vW,{t)) J] 11 + ^^') 

j=i j=i fc=i j=i 

where fii — — A^^'^-^ EU^'^ , i < d. For an analogical multinomial models (to those 



that we used for one Poisson process) we can prove a similar result to Theorem \2.1\ 
The proof for this case can he done in a similar way to the proof of Theorem \2.1\ 
and hy using the same ideas. For simplicity we provide the proof only for the case 
with one Poisson process. 



Error Estimates for American Options 



5 



Remark 2.3. By using the Cauchy-Schwarz inequality and the Chebyshev's in- 
equality it follows that we can set M{n) — [cn^^^] for some constant c ([x] is the 
integer part of x). Thus in the n-step mutinomial model which is given by 
the number of growth rates is proportional to n^/^. IfU^^'' takes on a finite number 
of values, then we can construct the multinomial models with a fixed number of 
growth rates by letting = for any i,n. In this case the proof of Theorem 

\2.1\ is simpler than for the general case and does not require Lemma [KM 

3. Auxiliary lemmas 

We start with a standard result, but since we could not find a direct reference 
its proof is given here for readers' convenience. 

Lemma 3.1. Let n € N and consider a probability space together with a filtration 
{^fe}fc=o '^'^'^ ^ positive d-dimensional adapted stochastic process {Q*-'^'}fc^o- 
any k let Y{k) = Q^*-^)) where (j)k : x ...W^ M+. Assume that 

E inaxoKkKn Y{k) < oo and define 

(3.1) A = esssupreME(Y{T)\gQ) 

where M is the set of stopping times with respect to the above filtration with values 

in {0,1, Assume that the sequence {( (k-i) , (fc^-i) )} "is a sequence 

Qi Qd k=i 

of i.i.d. random vectors such that for any k, (— (^rrTy, ^^-i) ) independent of 

Qi Qd 

Gk-1- Let p{ ) be the probability density function of {-^^ ^ -^^) . Then 

Qi Qd 

(3.2) A = Vo(g(°)) 

where -00 is given by the following dynamical programming relations 

(3.3) = 0„ and for <k<n, V'fc(a;*°\ xt^), x^) = max(</.fe(x("' , x^^) , 

^^''^)' X^eR" V'fc+i(a;*°\ a;(i\ x^^), a;'='^)p(zi, z„)dzi...dz„) 

where x^'^ = (xj'^, x^'^) € M!^ is given by x'l'^ = x'^^^ Zi, i < d. 

Proof. It is well known (see [13]) that A = A(0) can be calculated by the following 
dynamical programming relations 

(3.4) A{n) = Y{n) and for < fc < n, A{k) = max(r(fc), E{A{k + l)\gk))- 
By using backward induction, p.3p and p.4p we obtain that for any < fc < n 

(3.5) ^(fc) = V^fe(g("),QW,...,QW) 

and the result follows. □ 

Next, we derive several estimates that will be used in this Section. For any n G N 
and 1 < i < d denote, 

(3.6) S^{t) = Sm exp((r + /i. - E ■=! 4/2)^ + EU 

ji'^ ^ maxo<fc<„_i supkT/n<t<{k+i)T/n ISf i^) - SY{kT/n)\, 
= ELiSuPo<t<T \SY{i)\ and D = ELi suPo<t<T \Sk{t)\. 
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By using the inequality | exp(a;) — exp(?/)| < |a; — y| exp(max(a;, y)) we obtain that 
for any I < i < d 



(3.7) Jk'^ <D»'(\r + - ^ ti 4/21^^/ 



/n - 



E^=i kijl maxo<fe<„-i sup;,j,/„<t<(fc+i)j,/„ \Wj{t) - WjikT/n)\ 
Fix I < j < d. From the scahng property of the Brownian motion it follows 
i?f max sup \W,{t)~Wj{kT/n)A <nE{ sup \W,{t)\^)<^ 

\0<'=<"-l fcT/ri<t<(fc+l)T/n / 0<t<T/n ^ 

for some constant ci. This together with ^3.7\i and the Holder inequality gives 

(3.8) Ej24P<C2n-'/' 

i=l 

for some constant C2. Let T^ "- C T he the set of stopping times with values in 
{0,i,f ,...,T}. Set, 

(3.9) sup E{eM-rT)F{S{T),T)). 

Lemma 3.2. There exists a constant Ci such that for any n 

(3.10) 0<V- yji) < Cin-i/^ 

Proof. The inequality < V ~ Vn^"^ is obvious. Thus it remains to prove that 
V - Vj^^ < Cin-^/^. Fix n e N and choose e > 0. There exist t G T such that 

(3.11) y < e + £:(exp(-rr)i^(S'(r),r)). 

Define the random variable a = min{< e {0, ^, T}| t > r}. Observe that 

a e r^'" and T < cr < T + ^. Thus from it follows 

(3.12) V - yi^^ < e + E{exp{^rT)F{S{T),T)) - E{exp{-ra)F{S{(j), a)) 

< e + E\F{S{t), t) - exp{-rT/n)F{S{a), 
< e + E\F{S{T),r) - F{S{a),a)\ + ^-^E{F{S{a), a))<e+^il + ED) 
+LE{Eti \S^{r) - S,{a)\) + '-^{F{S{0), 0) + LD + LT{1 + D)). 
Next, set the event Q = {N{cf) > N{t)}. Notice that for any 1 < i < d 

N{t) 

(3.13) |5.(r) - 5,(a)| < IqD + 2( sup J] (1 + ul'^))J^^- 

0<t<T , 
J = l 

^From p.Sp . p.l2p - p.l3p and the Cauchy-Schwarz inequality we obtain that there 
exist constants C3 , C4 such that 

(3.14) k'^^ < e + C3n-i/4 + c4/P(Q). 

^From the strong Markov property of the Poisson process (with respect to the 
natural filtration generated by S) and the inequality t < a < t + we obtain 

AT 

(3.15) P{Q) < P{N{T/n) > 0) = 1 - exp(-AT/n) < — . 

n 

The result follows by combining p.l4p - p.l5p and letting e J, 0. □ 
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For any n define S^-'' = (5^", 5^'") by 
(3.16) 

d d N{t) 

Sl^-it) = SM exp((r + 4/2)^ + E '^^^^^i*)) IT + ^r), 1 < ^ < d. 

j=i j=i j=i 

Let 7"^'" be the set of stopping times with respect to the filtration {a{S{u), S^'"{u)\u < 
i)}?Lo with values in {0, ^, ^, ...,T}. Set 

(3.17) = sup E{exp{-rT)F{S^'"-{T),T)). 

Lemma 3.3. There exists a constant C2 such that for any n 

(3.18) \V^'^-Vj^'^\<C2n-'/^. 

Proof. ^From Lemma [3. II it follows that vi^-* = sup^-g-y-g.^ E{exp{—rT)F{S{T),T)). 
This together with (12. 5p gives 

(3.19) \V^^^ - Vl^^l < LE{supo<t<TEU - < L x 

LEiD^)Eij:u BuPo<*<T I nfJ(i + ui'^) nfJ(i + t^r'-'')i)- 

Next, we estimate the last term from p.l9p . Fix 1 < i < d. For any n,j € N 
set A:^'^' = [/f 'J' - The sequences are independent of the Poisson 

process N. Thus 

(3.20) i^(sup„<,<^ I nS^(l + U^'^) - U^fil + U-^')\) = 

EZi miT) = fc)£;(maxi<™<fc i n;ii(i + c/?'^) - 
nr=i(i + + < EZi miT) = k) X 

£/ ( maXl<„j<fc = l X/l<(;i <q2<...qj<m 

ni<.<™.^{,„...,,,}(i + f/i'^^)) < EZi miT) = fc) X 

E(^J2j=l J2l<qi<q2<...qj<m 111=1 l^i I I\l<s<m,s(^{qi ,...,qj}i^ + = 

Er=i^(^m = fc)E™=iEr=i (7)(i?iAr'i)-'(i +i?t/,r^)™-^- = 

from (12. 9p and the inequality \x"^ — y™| < m\x — ?/|(max(x, y))™^^ (for x, y > 0, 
and m S N) we obtain 

(3.21) = Er=i ^(^m = fc) E™=i((l + EU^'^ + E\A-''\r - 

(1 + EU^'Y) < n-'/' EZi PiNiT) - k) ELi ™(2 + >)— ^ 

= ri"^/^ EZ=i mP{N{T) > m)(2 + EU^^^)"'~^ < 
«"'/'E™=i "^^^(2 + SC/f ')"-i = ATexp(Ar(2 + ^C/f ^))n-i/8. 
By combining ((3l9| - ([3:2T|) we obtain ([3l^ . □ 
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For any nGN, 0<fc<n and 1 < i < d set 

(3.22) 5!*^ = niaxo<„<„ 1(1 + (1 - cxp{- XT /n))EUl''Y"' " exp(Ai,mT/n)|, 

2;n,k _ J2j^i^{N{jT/7i)~N{{j-l)T/n)>l} and 

S^--{t) = 5.(0) exp((r - E -=1 4/2)^^/^ + E -=1 ex., W^,(fcT/n)) x 

_ — uf:;'ii+ur) fo,M:<t<m)i:. 

(l + (l-oxp(-AT/n))_E(7r'') " " 

Notice that Z''^'' - N"^'' (the last term was defined before ^J^). iFrom ([O)) 
and (|2.9p it foUows that there exists a constant C5 such that 

(3.23) 5^:' < ^, VI < z < d. 

Denote S*^'" = (S*^'", 5^'")- Let T^'" be the set of stopping times with respect to 
the fihration {a{S{u), S^'^'iu), S^''^{u)\u < t)}f^Q with values in {0, |, ^, T}. 
Define 

(3.24) y„(3) ^ £:(exp(-rr)F(S'2^"(r),T)). 

Lemma 3.4. There exists a constant C3 such that for any n 

(3.25) |V;('^ -KP^I < C3n-i/«. 

Proo/. Fix n e N. Introduce the events Qk = {N{kT/n) - N{{k- l)T/n) > 1}, 
k < n. Notice that for any 1 < i < d 

(3.26) maxo<fc<„ jS'f "(/cT/n) - Sl'^'ikT/n)] < Iju^^^^Q,} supo<t<T(S',^'"(t) + 

5,2^"(i)) + G« exp(|A.,r|) supo<,<^ Sl-'-{t). 

^From Lemma [3.11 it follows that l^p-* = sup^g-T-a.n -E(exp(— rr)F(S'(r), r)). It is 
easy to verify that the terms i?(supQ<«j^(S'j^'"(i) + Sf'"'{t))'^) are uniformly (with 
respect to n) bounded. Thus (|3.26l) and the Cauchy-Schwarz inequality gives that 
there exist constants cg , C7 such that 

(3.27) \VS''^ - Fi'^l < LE j:ti niaxo<fe<„ \Sl'\kT/n) - 5^"(fcT/n)| < 

C6n-i/« + crVWk^iQk) < cen-'^^ + c,^^ 

and the result follows. □ 

Next, we derive estimates on a discrete probability spaces. For any n and < 
k < n define 

(3.28^f '"(t) = 5.(0) exp(fcrT/n + E™=i E ■=!( - 4^/(2"))) 



-X, for kT/n < t < {k + l)T/r 



[l + {l-cxp{-XT/n))E„ 

and yi^^ = sup^en En{cxp{-TT/n)F{T, ^^'"(r)), where 5^'" = (5^", 5^'"). 
Lemma 3.5. There exists a constant C4 such that for any n 
(3.29) \Vj^'^ ^Vn\<C\n-'/^ 
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Proof. For any 7i e N, < to < n and 1 < i < d set, 

(3.30) Br^" - ln(l + E -=1 ^/TJ^'^^,it^) ~ ^^=^{^'^^4^^ - ^hT/{2n)\ 

Observe that for any I < j < d, ICj^"*! < Vd+ I a.s. Thus by using the Taylor 
series of ln(l + x) we obtain 

(3.31) |B™'"-Cr'"| <C8n-3/2 

for some constant cg. ^From (I2.1ip it follows that for any i, n the random variables 
|^m,n|n^^^ ^^^^ ^ ^ with mean 0. This together with the Doob-Kolmogorov 
inequality (see [9]) gives 

m n 

(3.32) i?„( max | ^ < 4i?„(| ^ Cf-^ < 4ni?„(|C,^-n < ^ 

l<m<n ^ — ' ^ — ^ 71 
^ ^ /c=l fc=l 

for some constant cg. It is easy to verify that the terms En{J2i=i suPo<t<T('5'i^'"(*) + 
S'j"^(t))2), n e N are uniformly bounded. From (P3T|) - (P3^ and the Cauchy- 
Schwarz inequality we obtain that there exists a constant C4 such that 

\vJi'^ ~ Kl < LE^iEti BuPo<t<T - 5(")(t)|) < 

i£;„(Eti niaxi<„<„ I ELi ^r'"! supo<,<^(5f '"(t) + st\t))) < 
LE^iHUiom-"^ + maxi<™<„ | ELi ^"^"1) supo<,<j^(5f '"(i) + S^^^ it))) 

□ 

4. Proof of main results 

The following result which we state without proof was established in [M] under 
more general assumptions (see Theorem 2.1, Corollary 3.1 and Lemma 5.1 there). 
This result is the main tool that we use in order to complete the proof of Theorem 
EH 

Theorem 4.1. Consider a probability space together with a sequence 

of i.i.d. d-dimensional random vectors and a d-dimensional random vector Y . 

Assume that EJ2Lii\xl^^? + \Y^\^) < 00, EX'^^^ = EY and 

(4.1) E{X^'^X^/~>) = E{Y^,) e {!,..., d}. 

For any z > it is possible to extend our probability space to (il, F, P) which 
contain a sequence of i.i.d. random vectors Y'^^\Y'^'^\ ... such that Y^^"^ ^ Y , and 
for any n G N 

(4.2) P( max ^ | ^ xt^ Y^\ > z) < ^Y.^{\x\'\' + \Y,\^) 

~ ~ i—1 ni — 1 i—1 

for some constant C which independent of X^^\ X^-^^ Y and z (E denotes the 
expectation with respect to P). Furthermore, for any fc > 1 the random vectors 

xw,... 

^ j^(;e-i)^y(fe)^y(fe+i)^ _^ are independent. 
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^From Lemmas I3.2H3.5I it follows that there exists a constant C5 such that for 
any n & N, \V — Vn\ < C^n^^^^ + \vJi'^^ — |- Thus in order to complete the proof 
of Theorem 12.11 it remains to establish the following lemma. 



Lemma 4.2. For any e > there exists a constant such that for any n 

(4.3) \VP ~V}^^\<C,n^-'^. 

Proof. Fix rt G N and e > 0. We start with proving the inequality 

(4.4) v;p) - i^i^) < Cen^-i 

^From ()2.1ip and Theorem l4.1l it follows that we can construct on the same probabil- 
ity space (51, F, P) two sequence of i.i.d. random vectors X*^^' , X^"^ F^^^ F^") 



such that - W{1), - V i^^'^ and 

(4.5) \^t^Jl\ E - > ^'"") ^ C^^'"" 

i—1 m— 1 

for some constant C^. Furthermore, for any k > \ the random vectors X'^^\ X^^~'^\Y'^^^ 

are independent. We can extend the constructed probability space such 
that it will contain also two independent sequences of i.i.d. random vectors F'^-* , i^^"^ 
and G(i),...,G("' which are independent of X^^), F^^), and satisfy 
- G(i) - (/9"'\ [/"'I were defined in Section 2). 

For any p < n and i < d set, 

(4.6) Z(P) = Y.U F^^\ A = maxo<,.<„ 

= 5,(0) exp(rr) maxo<fc<„ e^TL^i E -=1 ('^^j-^j"^ " ^ljT/{2n))) 

and Dj = SM exp(rr) maxo<fe<„ exp(E^^, Y^Ui^^^J^o"'^ " 4^/(2")))- 

It is easy to verify that there exists a constant C and for any p > 1 there exists a 
constant Cp (the above constants does not depend on n) such that 

d d 

(4.7) ^ < G and i?((^ Af + i^f)'') < G^. 

i=l i=l 

Define ^(t) = (^i(t), ^^(i)) and S{t) = (^i(t), ^^(t)) by 

(4.8) S,{t) = 5.(0) exp(rfcr/n + YL=i E ■=i('^^^-^r' " ^lT/{2n))) x 
^^'=^ ^ + - and S'.(i) = S•.(0)exp(r/^T/r^ + 



(l+(l-cxp(-AT/n))£'GP'y 



nfi'i'(i+G^^') 



^+(l-exp(-AT/ra))£;Gf' ' 



where l<i<d, 0<yfc<nand^<i< IMllI^. 

Next, let T be the set of stopping times with respect to the filtration {a{Si{u)\u < t)}^^^ 

T TT 
n ' n 



with values in {0, T} and let T be the set of stopping times with respect 
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to the filtration {(7{S{u), S{u)\u < i)}t=o ^^i^^ values in {0, ^, ...,T}. Set, 

(4.9K - sup^^f E{eM-rr)F{SiT),T)), V = sup^^^ EieM~rT)FiS{T),T)) 

and Q = {maxi<,<„EtilELi^l"->^/'"^l >-~^/n. 

Observe that S - 5*2'" and 5 - S^^". ^From Lemma O it follows that V = vi^^ 
and f = Vji^\ in the last equality we used the fact that X^^\ X^'"'^^ , F^^) , F^") 
are independent for any A: > 1. Since T C T then from (|4.5p . (I4.7I) - (|4.9I) . the in- 
equality I exp(x) — exp(y)| < |a; — y| max(exp(a;), exp(y)) and the Holder inequality 
we obtain 

(4.10) - Vji^^ < sup^^f EieM-rT)FiS{T),T) - 

sup^^f E{exp{^rT)F{S{T),r)) < i^(supo<,<^ Eti \S^{t) - km 
< LE{{n-^'^ niaxi<,,,<d + lQ)(Eti D^{Df + Dj))) 
< i(7((C6)i/fn-i/8+e(C',)i/9 + maxi<,,,<d \a,j\Cin-^/^) 

where p — and q — This completes the proof of the inequality (|4.4I) . 

The inequality Vn — Vn < C^n'^~'s can be proved in a similar way, just take 
^ and yd) ^ ly(^). □ 

Remark 4.3. In several cases our approach can be extended also for path depen- 
dent payoffs which satisfy Lipschitz type conditions. Let M[Q,t] be the space of 
Borel measurable functions v = {vi, ...,Vd) ■ [0,t] M'' with the uniform metric 
dot{v,v) = supo<,<tE^=i |w»(s) 

—Vi{s)\. For each t > let Ft be a nonnegative function on M[0,t] such that for 
any t > s > and v,v d M[0, t], 

\F,{v)~Fs{v)\<Ldos{v,d) and \Ft{v) - Fs{v)\ 

< L{\t - s|(l + sup„g[o.t] Eti + sup„g[,^t] E.f=i - Vi{s)\). 

Consider an American option with the payoff process Y(t) ~ Ft{S) where S = 
^(a;) G Af[0,oo) is a random function taking the value S{uj,t) at t d [0,oo). 
When considering Ft{S^) for t < oo we take the restriction of S to the interval 
[0,t]. The term V = sup^g-y- i?(exp(— rT)y(r)) gives an arbitrage- free price for our 
model. For the multinomial models we consider the payoffs F*^")(fc) = Fkx (5*^"-*) 

and the arbitrage- free prices ¥„ = sup^^^^ £'„(exp(— TT/n)y'^")(T)). Lemmas \3.2l 
\3.5\ and \4.2\ can be extended to this setup in a way that does not ruin the esti- 
mates of Theorem \2.1l The problem is with Lemmas 1 3. 3[\S. 4\ For path dependent 
options the equality before l\3.19^ does not follows from Lemma \3.2\ ( although is 
correct) and the the first inequality in \3.2T^ should be replaced by jvi"^^ — < 
i£' EiLi suPq<j<7- \Sf'"{t) — Sl'"'{t)\. The way to fix it is to consider a piecewise 
constant approximations of S, S^"^{t) = Skx/n for kT/n < t < {k -\- l)T/n, 
< k < n. If we could provide an estimates of the term i?niaxo</j<„ \FkT (S) — 

FkT {S^'^^)\ then both of the Lemmas \3.S\^3^\ could be extended. Observe that, 

S(maxo<fe<„ \F!^{S) - < iS(supo<i<T Eti \S^{t) - S^r\i)\)- For 

the Black-Scholes model (no Poisson process) the last term was estimated in I13.8\) . 
Thus for the Black-Scholes model the estimates from Theorem \2.1\ remain valid for 
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path dependent options which satisfy Lipschitz type conditions. However, if we allow 
jumps of compound Poisson type than in general the term £'(supo<f<y X^^Li 

S\"'\t)\) should not tend to 0. Consider a specific type of path dependent options 
which are given by Ft{S) = niax(M, suPq<;«2^ niaxi<i<£; where M > is 

some constant (Russian options). Recall the terms Jn^ that were defined in h3.0\) 
and the events Qk which were defined before h3.26\) . For Russian options we have 
the following inequality 

niaxo<fc<„ \F^{S) - < iI{uj"^,Q,} supo<,<j, f^t(5) + 

BuPo<.<TnS^(l + C/i^'^)Eti^^^'- 
By using estimates that were derived in Section 3 we obtain that -E(inaxo<fc<„ \FkT_ {S) — 
FkT (S**-"-*)]) is of order n~^/'^ and so the estimates from Theorem \2.1\ are also valid 
for Russian options in Merton's model. 
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